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RESEARCH MEMORANDUM 

PRELIMINARY MEASUREMENTS OF TURBULENCE AND TEMPERATURE 
FLUCTUATIONS BEHIND A HEATED GRID 
By A. L. Kistler, V. O'Brien, and S. Corrsin 

SUMMARY 


Preliminary measurements have been made of velocity and temperature 
fluctuations in the flow behind a heated grid in a uniform airstream. 
Temperature correlation shows a reasonable degree of isotropy, and the 
temperature fluctuations die out at large distances more slowly than the 
turbulence, as has been predicted theoretically under some strongly sim- 
plifying postulates. 


INTRODUCTION 


A comparison between the correlation equations for concomitant 
(incompressible isotropic velocity and temperature fields showed marked 
differences attributable to the fact that velocity is a vector while 

temperature is a scalar (refs. 1, 2, and 3) . 1 The relative decay rates 
were calculated voider strongly simplifying assumptions , and the rela- 
tive "microscales" and relative wave numbers characterizing the very 
fine structures were deduced in terms of the Prandtl number of the fluid. 
Since the theoretical predictions are essentially conjectural, it is 
necessary to get at some facts through measurement. 

A hot grid with thermal mesh equal to momentum mesh is as simple 
as any, so this configuration was selected, with the expectation that 
approximate equality of integral scales would be obtained. 

The analysis of reference 1 suggests that the flow with equal inte- 
gral scales for velocity and temperature fields may not be the simplest 
case analytically, but could be the si mplest for experimental realizatio n. 

-*-As pointed out in references 1 and 2, the work applies equally well 
to isotropic turbulent mixing between two different components provided 
that the molecular mass transfer coefficient is nearly constant over the 
concentration range encountered. 
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This report describes the preliminary data. It is hoped that with 
improved grid and measuring procedure, less scattered results will be 
gotten in the measurements now in progress . 

The investigation was conducted at the Johns Hopkins University 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 

SYMBOLS 


% (*l) 

one -dimensional power spectrum of u-fluctuations 

E -a ( k i) 

one -dimensional power spectrum of -8 -fluctuations 

f 

"longitudinal" double velocity correlation, 
u(x,y,z)u(x + |,y,z) 
u' (x)u' (x +| ) 

g 

"lateral" double velocity correlation, + 

u t2 (x) 

k i 

wave number in x-direction 

k c,3 k c 

spectral "cut off" for Ep ^k and E^ , respectively 


integral scales of g, f, and m, respectively 


grid mesh size for momentum and heat =1 in.^ 

m 

temperature correlation 

r 

space interval 

t 

time 

u 

turbulent velocity fluctuation in x-direction 


Cartesian space coordinates; x is alined with mean flow 

7 

thermal dlffusivity coefficient , K/pCp, where K is 
thermal conductivity coefficient, p is density, and 
Cp is specific heat at constant pressure 
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n 

•9 

A, Ag 


V 

I 

a 


space interval in y-direction, =Ay 
temperature fluctuation 

dissipation scales for velocity and temperature fluctua- 
ti ons , respectively 

kinematic viscosity coefficient 

space interval in x-direction, =Ax 

Prandtl number, v/7 

EQUIPMENT AND PROCEDURES 


Approximately isotropic turbulence is produced at some distance 

behind a square -mesh biplane of grid of ~ - inch round wooden dowels 

4 

spaced 1 inch on centers. The grid is set at the upstream end of the 
working section of the wind tunnel sketched in figure 1. 

As a first method of generating isotropic temperature fluctuations 
in this turbulence, a grid of Nichrome heating wires was introduced about 

li inches upstream of the dowel grid, with wires geometrically in phase, 

that is, alined so that their wakes passed over the rods (fig. 2). In 
this way the hot air passed through the turbulence -generating shear zones 
of the main grid, so that temperature and velocity fluctuations were pro- 
duced in the same local neighborhoods. 

Turbulence and temperature fluctuation measurements were made with 
the hot-wire anemometry equipment described in reference k. The experi- 
mental procedures for simultaneous velocity and temperature fluctuations 
are given in reference 5* 

The measurements were all made at a mean speed of 1^ feet per sec- 
ond and a mean temperature rise (across the grid) of about 5° C. The 
actual heating wire temperatures were about ^>00° C. The resulting tem- 
perature fluctuations were inconveniently small, but this insured neg- 
ligible influence of density variations upon the fully developed turbu- 
lent velocity field - as attested by both turbulence level and correlation 
measurements . 

Because of the extremely small temperature fluctuations it was nec- 
essary to use unusually high resistance hot-wires, 0.0001-inch platinum 
of about 90 chms resistance. This wire was set in the form of a V (with 
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plane perpendicular to mean flow) in order to reduce the error due to 
finite wire length. The high resistance necessitated correction for 
heating current fluctuations when the system was operated as an anemometer. 

As has been found before (ref. 6 ), the velocity field was not quite 
isotropic within the precision of measurement: v f « 0 . 9 u f at x/m = 100 ; 

but this is apparently insufficient anisotropy to cause difficulties . A 
reasonable degree of isotropy in the thermal field is indicated by the 
approximate equality of longitudinal and lateral temperature correlation 
functions ( fig . 5 ) • 

For convenience the spectrum was measured without heating, after it 
had been determined that the u' -levels and correlations were not appre- 
ciably changed by the heating. 

The temperature spectrum was measured with "mixed” sensitivity, that 
is, wire temperature set for the same order of response to velocity and 
temperature. This was done because under resistance -thermometer condi- 
tions (heating current < 1 milliampere) the signal was of the order of 
the amplifier noise level over most of the frequency range. 

With mixed sensitivity the hot-wire responds to both velocity and 
temperature spectra, which must be "separated." Simple superposition of 
the two energy spectra follows from the assumption that all "harmonics" 
of the two spectra are uncorrelated. Since u^fl =0 in this field, such 
an assumption seems reasonable. Then the temperature spectrum is obtained 
by subtracting the (unheated grid) velocity spectrum from the mixed 
spectrum. 

A Hewlett-Packard model 300 A wave analyzer was used as variable 
frequency band pass filter. 

No correction has been made for wire length. 


EXPERIMENTAL RESULTS 


The velocity and temperature decay curves are plotted in figures 3 
and 4, respectively. A number of runs are given for each to indicate 
the repeatability of the data as well as the equality of hot- and cold- 
grid flows (fig. 3 ) and the consistency of two difference measuring 
methods (fig. 4). The upvalues are a bit lower than those obtained 
behind a similar grid in a different wind tunnel (ref. 6 ). 

Figure 5 shows longitudinal and lateral temperature correlations 
at x/M = 69 . Of course it was necessary to keep the downstream probe 
out of the wake of the upstream one, so the Ax coordinate is actually 
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inclined to the x-axis by about 5°- The two sets of points match fairly 
well, indicating that the temperature field is not far from isotropy. 

Figures 6 and 7 compare these curves with the corresponding velocity 
correlations. The close agreement of the longitudinal cases is especially 
noteworthy . 

The experimental scatter is too great to permit a determination of 
the dissipative scales ("microscales") by graphical fitting of parabolas 
to the vertices of these correlation data. However, the integral scales 
have been calculated: 





>(i) 


where f(r) and g(r) are the Karman-Howarth notation so that 


_ u (*,y, z)u{x + g,y,z) 
u' (x)u' (x + I ) 


g ( n ) _ u(x,y,z)u(x,y + t),z) 

Also 


s ' | 9(x,y,z)^(x + g,y,z) 
•3 ' (x)3 ' (x + | ) 
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and 


di(ti) = + ti,z) 

* 2 (x) 


are the longitudinal and lateral temperature correlation coefficient 
functions . 

The results are 


L = O.kO in. 


= 0.70 in. 




0.66 in., 
0.68 in . , 


laterally 

longitudinally 


> ( 2 ) 


In passing one notes that — = 0.57 i- O.5O, the value corresponding to 

L 1 

exactly isotropic turbulence. Since M, the momentum grid mesh, and 
M$, the heat grid mesh, are both 1.0 inch, the values in equations (2) 
are also l/m, Lq/M, and L0/ML0, respectively. 

The time spectra of *8 and u at x/M = 23 were determined as 
described in the previous section and are given in figure 8. These may 
be considered good approximations to the one -dimensional longitudinal 
space spectra, except perhaps for the lowest wave number range (refs. 7 
and 8). The sketched curves indicate a possible trend, but should not 
be taken seriously; within the experimental precision there is no sig- 
nificant difference. The circuit noise has been subtracted out for both 
sets of data. 

These spectra have been used for approximate calculation of the 
dissipative scales of velocity and temperature fields: 



JL. 

X 2 


(3) 
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(4) 


The results of graphical computation are 


A^ = 0.26 in. 

A = 0.19 in. 




(5) 


at x/M = 23. 

With the assumption that velocity and temperature fields are suf- 
ficiently near to isotropy A(x) and A^(x) can be calculated from 

the intensity traverses, as indicated in the next section. The resulting 
faired curves are plotted in figure 9> with the above values spotted in. 


COMPARISON WITH THEORY 


Of course the theories of isotropic turbulence and of isotropic 
scalar fluctuations in isotropic turbulence are still unsolved; in fact 
they have not even been formulated in a determinate way. Nevertheless, 
a number of consequences of isotropy have been inferred, including some 
particular (and relatively simple) forms of the averaged differential 
equations. Of especial interest here are the so-called ’’decay equations.” 
For isotropic fields decaying in time, there are (refs. 1 and 9) 



( 6 ) 


arid 


aa 2 

at 



(7) 


The actual experiments are done in steady state behind a grid with 
the various averaged quantities depending upon x instead of t. Taylor 
(ref. 9) has pointed out, in effect, that if these vary sufficiently 
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slowly with x equation (6) can be applied to the grid flow with the 
substitution x/U = t. 2 The same holds for equation (7), so that 


du 2 

dx 


-10 * 

u x 2 


d-9 2 

dx 



( 8 ) 

(9) 


These two equations permit calculation of X(x) and X^(x) from 

the experimental data on u'(x) and d'(x), as mentioned at the end of 
the previous section. The faired results are drawn in figure 9. The 
appreciable disagreement at x/m = 23 between X's calculated from 
decay and from spectra is indicative of inaccurate spectra (the less 
accurate of the two kinds of data) and/or of appreciable deviation from 
isotropy this close to the grid. 


The ratio X3/X varies almost linearly from 1.23 at x/m = 25 

to I.58 at x/M = 100 (fig. 9)* The only other measurement of this 
ratio was on the axis of a heated round turbulent jet (ref. 4) where 
the value was about 1.34. * 


In reference 1 it is shown that for two special "kinds" of isotropic 
fields this ratio has the same simple value: 



( 10 ) 


where a is the Prandtl number v/7. The two kinds of isotropic fields 
are: 


(a) Reynolds and Peclet numbers are so small that the convective 
effects are negligible for both heat and momentum. Therefore this is 
not. really turbulence in the ordinary sense . 

(b) Reynolds and Peclet numbers are both very large and each field 
is assumed to have a completely "self -preserving" correlation function 
during decay. In this model, due to Von Karman (ref. 10), the dissipation 

2 ' ~ — 

In fact, ^this requires perhaps more justification than given by 

Taylor or by Karman and Howarth. This has been discussed at a seminar 
lecture at M.I.T., May 1952, and printed copies are available from the 
Aeronautics Department of The Johns Hopkins University: "Remarks on the 

Comparison of Grid-Produced Turbulence With the Theory of Isotropic 
Turbulence, " by S. Corrsin. 
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and conduction terms of the two correlation equations are neglected, but 
decay is indirectly included by substitution of the decay equations into 
terms of the correlation equation. 

Equation (10) can also be deduced by an approach somewhat different 
from the above: 


Assume the Reynolds and Peclet numbers so high that the ’’inertial 
subranges” ^with E^k^ and E^k-^ both ~ / are sufficiently 

extensive to permit the approximations (ref. 2) 


N 


-5/5 


when < k n < k 


o 1 


= 0 elsewhere 


( 11 ) 


( k l) ~ k l" 5/ ^ when •3 k o < k ! < * k c 


= 0 elsewhere 




The ’’cut off” wave number k c is just the inverse of Kolmogoroff f s 
microscale. 


k c = — 
c n 


( 12 ) 


and. ^k c is the analogous thermal quantity (ref. 2) 


k _ 2n 
* c ■ 


( 15 ) 


In reference 2 it was shown that 


i3 k c = 0 5/4 

k c 


(14) 


Substituting equations (ll), (13)> and (l^-) into ( 3 ) and (4) and 

-relative to k_^/3 and ^k 0 ^/3 relative to ^k c ^/3^ 


neglecting k Q 1 J relative to k Q 
there results 


and 
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( 15 ) 


If one now restricts to fields in which the thermal large structure is 
roughly equal to the velocity large structure (as in these experiments), 
k Q ~ and one arrives again at equation ( 10 ). 


In fact, none of the three sets of assumptions leading to equa- 
tion (10) is directly applicable to these experiments. The Reynolds 
and Peclet numbers are not low enough for case (a) and not high enough 
for the other two cases. Furthermore, it is known that the correlation 
functions of velocity do not remain completely similar during decay 
(ref. 11). The important conjecture to be made is that a result encoun- 
tered under such divergent sets of assumptions may have approximate 
validity in a wide range of situations - perhaps including the actual 
one . 


Taking o = O .72 for air, equation (10) gives 

T “ 1 ' fiT 

which is in reasonable agreement with the downstream values computed 
from the decay data (fig. 9 )- 

The A^/A curve is just a particular way of displaying the decay 

data. Two other ways are instructive: Figure 10 shows the variation 

in Reynolds and Peclet numbers 




> (16) 


The apparent constancy of P^, is evident. Figure 11 shows the relative 

decay rates of temperature fluctuations and velocity fluctuations . For 
small values of x/m they die out at about the same rate, but for larger 
values of x/M -9 ' dies out at a decidedly lower rate than does u' . 

For fully isotropic fields the ratio of equation ( 7 ) to equation (6) gives 


R - U ' A 

r a = — 


% E 


U'A, 
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( 17 ) 


and if equation (10) holds 


which is very close to the measured results at large values of x/M. 

This is just another way of showing that A^/A does approach the value 

in equation (10). 

In the absence of more complete data, no attempt has been made to 
compare the individual decay rates of i3 1 and u' with the theoretical 
predictions that can be obtained by a variety of postulates. 


The principal conclusions to be drawn from the foregoing results on 
measurements of turbulence and temperature fluctuations behind a heated 
grid are that: 

1. It is experimentally feasible to set up and measure concomitant 
velocity and temperature fluctuation fields downstream of a heated grid. 

2. The two fields are moderately isotropic sufficiently far 
downstream. 

3. Also, sufficiently far downstream, the dissipation scale ratio 
is approximately given by 


a value deduced theoretically under widely divergent postulates, none of 
which corresponds closely to the present experimental conditions. 


CONCLUSIONS 
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4. As a consequence of conclusions 2 and 3 , far enough from the 
grid the temperature fluctuations die out more slowly than velocity 
fluctuations : 


4’(*)~ [u'(xj] 5/5 

corresponding to equation (l8) . 

The Johns Hopkins University, 

Baltimore 18, Md., January 13, 195^- 
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Figure 1.- Sketch of wind tunnel. 
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Figure 2.- Hot grid photographed from downstream face. Parts of square 
heating mesh of double wires can be seen through wooden dowel grid. 
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Figure 3*- Turbulence decay. 
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Figure k.- Temperature fluctuation decay. 
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Figure 5*- Longitudinal and transversal temperature correlations. 
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Figure 6.- Transversal velocity and temperature correlations, 
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Figure 7.- Longitudinal velocity and temperature correlations. 
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Figure 9*- Vs calculated from u 1 and $ ? decay curves (points from 

measured spectra) . 
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Figure 10.- Reynolds and Peclet numbers. 
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